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(~| . Abstract. The characteristic polynomial for a chain of dipole-dipole coupled 

Oh| two-level atoms with nearest-neighbour and next-nearest-neighbour interactions is 

developed for the study of eigenvalues and eigenvectors for single-photon excitations. 
Such a system is mathematically equivalent to an XX spin chain in an external 
^ ' magnetic field. We find the exact form of the polynomial in terms of the Chebyshev 

O^. polynomials of the second kind that is valid for an arbitrary number of atoms and 

coupling strengths. We then propose a technique for expressing the roots of the 
04 ' polynomial as a power series in the coupling constants. The general properties of 

J^ , the solutions are also explored, to shed some light on the general properties that the 

QQ I exact, analytic form of the energy eigenvalues should have. A method for deriving the 

0^ ' eigenvectors of the Hamiltonian is also outlined. 
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Vh ' 1- Introduction 



The properties and dynamics of linear chains of atoms under the nearest-neighbour 
couphng approximation have been the subject of considerable interest for many 
decades [H [21 [3] . Recently, a new interest has arisen in linear chains in the context of 
quantum information and computation [Ij, where linear chains have been proposed as 
entanglement carriers that can be relatively easily realized in practice by diverse physical 
systems such as cold atoms trapped in optical potentials [3, E], quantum dots [3, El E] 
and Josephson junction arrays [lOj. Despite the simplification that takes into account 
only the nearest neighbour interaction between the atoms, the analysis of the properties 
of a linear chain containing a large number of atoms are still complex. Therefore most 
of the developments have been focused on the analysis of the energy and entanglement 
of the ground state only. Numerical results for the ground state energy and two-particle 
correlation functions were obtained using the modified Lanczos method [HI [12] . Some 
approximate analytic results for the ground state energy and higher excitations were 
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also reported [TB] . In addition, the thermodynamic properties of a frustrated Heisenberg 
chain were studied [H] and fluctuations in the axial Ising model were analysed [15]. More 
recently, numerical investigations were conducted [IE] , and also the scahng of excitations 
were analysed [T71 [T8] . The work of Lamas et al. [19] have introduced approximations 
for handling a spin chain with different coupling strengths along the chain. Apart from 
numerical results, an analytic approach of finding the ground state was proposed [20] , 
and the properties of frustrated Heisenberg chains were also studied |21j . 

As we have already mentioned, exact treatments of linear or closed chains are 
limited to the studies of energy and entanglement of the ground state or to a small 
number of atoms of an arbitrary excitation [22l [23] . In the present work, we consider 
the dipole-dipole chain, which is equivalent to an XX chain in a magnetic field, and 
derive the characteristic polynomial equation for the energy states and investigate the 
properties of its solutions. Derivation of the eigenstates of a chain of interacting atoms 
is of great importance for the study of entanglement creation, its transfer through the 
chain and stability against the decoherence. We also go beyond the nearest-neighbour 
coupling approximation to include the next nearest-neighbour coupling. 

It is generally believed that the nearest-neighbour approximation is a correct 
approach to study the properties and dynamics of linear chains. However, for realistic 
physical implementations of linear chains some directions of investigations are important 
such as the role of long range dipole-dipole interactions. It is the main purpose of 
this paper to address this issue by considering effects of the next-nearest-neighbour 
interactions. We shall deal with the eigenvalue equation for a chain of n coupled atoms 
and extract the characteristic polynomial for the eigenvalue equation of the Hamiltonian 
of the system for a single excitation. We then propose a technique of expressing the 
roots of the polynomial in terms of a power series in the coupling constants. Deriving 
the exact characteristic polynomial moves the research field closer to the goal of exact, 
analytic results for a chain of arbitrary length and arbitrary excitation, for it is the roots 
of this polynomial which give the energy eigenvalues of the chain. 

2. System Hamiltonian 

We consider a linear chain of n identical two-level atoms. The atoms are assumed to be 
equally spaced and confined to fixed positions, and each atom can only interact with its 
nearest neighbours through the dipole-dipole interaction. A schematic of the system is 
shown in Figure [U 

The dynamics of the system are governed by the Hamiltonian: 

n 

H = J2 ^oSl + a Y, S^Sr + b J^ S+S-, (1) 

where the parameters a and b describe the nearest-neighbour and next-nearest-neighbour 
coupling constants, respectively. The operators S^ and S^ raise and lower the state 
of the i-th atom by a single excitation, and S^ = S^S^^ . This definition of S^ makes 
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Figure 1. Physical layout of the linear chain. The atoms are equally separated from 
each other by a distance d, and are polarized in the same direction in their dipole 
moments /x and their two- level energy structure (shown beneath each atom). 



the energy of the ground state of the system equal to zero. We first briefly show that 
the above Haniiltonian is equivalent to an XX chain in a magnetic field [Tl[2j. We first 
express the raising and lowering operators in terms of Pauli matrices: 



St 



[(tI + ktI'' 



'k - 2V^fe -r '-^fc;, Sf, -^{(7^-icrl). (2) 

These relations identify the internal states of a two-level atom with the z component of 
spin for a spin-| particle. Substituting these relations into ([I]) and making use of the 
commutation relations for Pauli matrices gives the analogous Hamiltonian for a spin 
chain: 



1 ^ 



H = -nu, + -u, Y^ al + ^a J] (af aj + af aj) + ^& E (^^^^i + ^'^P- (3) 



Ignoring the first term, which just shifts the energy scale, the second term represents the 
coupling of an external magnetic field to the spins of the particles, while the second and 
third bracketed terms represent XX coupling for nearest and next-nearest neighbours. 

It is important to note that our work does not depend on specific features of the 
dipole-dipole interaction; all of our results are therefore also applicable to spin chains. 
We will, however, give specific examples in terms of the dipole-dipole coupled chain to 
fix our ideas, since it provides an easy way to think of a system described by couplings 
between raising and lowering operators. 

Although we will regard a and h as arbitrary parameters, their magnitudes can be 
explicitly determined, for the case of dipole-dipole coupling, from the expression p^[25] 






-r 



[1 






smxj 



+ 



cos Xi 



A~ \ ^~ \r- ' -J I i ^ • V- -\r- - ojj i ^2 ' ™S I ' V / 

where Xij = 27rrjj/Ao is the distance between the atoms i and j relative to the resonant 
wavelength Aq = c/2ttujq of the atomic transition, F is the decay rate of the atomic 
transition, assumed the same for all atoms, and fi and Vij are unit vectors in the 
direction of the atomic dipole moments and the displacement r^j between the i-th and 
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Figure 2. Dipole-dipole interaction strength as a function of ry/Ao for two different 
orientations of the atomic dipole moments relative to the interatomic axis, (il-fij (solid 
line) and (i \\ Vij (dashed line). 



j-th atoms. The dependence of Qij on the distance between two arbitrary atoms in the 
chain is illustrated in Figure El In the present case, a = fii,i+i and b = fii,i+2- 

We will first derive the eigenvalue equation for the Hamiltonian and investigate 
is its symmetry properties. This will be followed by a discussion on the form of the 
eigenstates of the Hamiltonian. 



3. Eigenvalue equation 

Since our system comprises of n two-level atoms, it is clear that up to n photons may 
be present in the chain. We will consider only one-photon states due to the complexity 
of the problem for higher-order excitations. Single-photon states may be useful in 
approximating higher-order states. 

Our goal is to find analytical expressions for the one-photon eigenstates and their 
energies, which are found by solving the eigenvalue equation H\iIj) = E\iIj), where H is 
given in ([T]). We proceed by writing the characteristic equation 



Tn = det{H-EI) 
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where X = ujq — E and / is the n x n identity matrix. The matrix T„, which represents 
the eigenvalue equation, is a n-dimensional five-diagonal matrix with two off-diagonals 
set by the nearest-neighbour coupling and the other two off-diagonal set by the next- 
nearest-neighbour coupling. For clarity, all entries in a determinant which are zero are 
simply omitted, and the order of the determinant (in this case, n) is placed as a subscript 
on the bottom right-hand corner of the determinant. 

The eigenvalues E are then given by the roots of T„, i.e. those values of E which 
make Tn{E) = 0. At this point, we could pick a specific numerical value for n and solve 
for the roots of T„ by expanding out the determinant. However, we would like to keep 
our results as general as possible; we therefore seek an alternate expression for T„ such 
that it is an analytic function of n. Such a form has the advantage that, as evidenced 
by the final result, one does not need to construct a matrix to find the eigenvalues. This 
can be an advantage particularly when n is large since less computer memory is required 
for the calculation. 

To find a form of T„ that is analytic in n, we first expand the determinant appearing 
in ([5]) to find a recurrence relation for T„. An expansion of the determinant leads to a 
closed recurrence equation 



T„ = AT„_i - aDn-i + haDn-2 - &^AT„_3 + 6^T„_4, 



(6) 
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(7) 



b a X 

The next step is to obtain an expression for T„ that does not have any D„ terms in 
it. We do this by finding a recurrence relation for D„. Having then two simultaneous 
recurrence relations involving T's and D's, we can then eliminate all the D's in the 
recurrence equation ([6]) leading to a recurrence relation solely in terms of T's. 

Thus, expanding the determinant Dn along the first column gives 

Dr, = aT„^i - bD^r^-i. (8) 

It is seen why we insisted on not expanding out the D's terms in ([6]). This procedure 
would be equivalent to repeatedly substituting ([8]), which would not end until all the 
rows of Tn have been Laplace expanded. This would give us an effectively useless n-term 
recurrence relation since the number of terms grows with n. 

We now use ([8]) to eliminate the -D„_i term from ([6]). This step gives 

Tn = XTn-l + 2baDn-2 ' o'T^-s - b^XTn-3 + b^Tn-4. (9) 
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Next, by adding 6T„_i to both sides of the above equation, we find 

T, + 6T„_i = A(r„_i + 6T„„2) + 26a(D„_2 + &^„-3) 

- a2(T„_2 + &T„_3) - &'A(T„_3 + 6T„_4) 

+ 6^(r„_4 + &T„„5). (10) 

Since -Dn-2 + bD^-z = aTn^^, as it is seen from ([8]), we finally arrive to the recurrence 
relation for T„: 

T„ = (A - 6)T„_i + (A6 - a2)T„„2 + ia% - b^\)T,,_s 

+ {b^-b^X)Tn-4 + b'T^.,. (11) 

The recurrence relation (ITT|) looks complicated and does not reveal any similarities to 
the well known recurrence relations for special functions. Therefore, we will try to find 
a non-recursive form for T„. 

4. Characteristic equation for Tn 

Since the recurrence relation flTTj) is a linear homogeneous equation with constant 
coefficients, we let T„ = x" to form the characteristic equation 

x^-{X- b)x^ - {Xb - a^)x^ - {a^b - b^X)x^ - (6^ - b^X)x - b^ = 0. (12) 

The roots of the characteristic equation can be easily found, and are of the form 

b 

x=< x± = lx+±Wxl-^b^ , (13) 



where 



and 



y±^lx-±Wx--^b^ 



X± = ^ (14) 



A = v/(A + 26)2 - 4a2. (15) 



Thus, with the roots flT^ . the general solution to the recurrence relation fITT]) can be 
written as 

T„ = G6" + P+< + P_x'!+Q+y'^ + g_y^, (16) 

where the parameters G, P±, and Q± are arbitrary constants. These constants are 
determined by requiring that flT6|) gives the correct result for five different values of 
n. We chose n = 1 through to 5 inclusive (the specific choices for n does not matter; 
the constants will not change). This gives us a set of five simultaneous equations or 
"boundary conditions" in the five unknowns G, P±, and Q±, which is sufficient to give 
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unique values for these parameters. These boundary conditions force the five arbitrary 
constants to take on the following values: 

-262 
^ "46A-862-a2' ^^^^ 

_ \{\ + 2b + A)-2a' X-, + b &[2&(A - 26 - 3x+) + a'] , . 
^ 2Ax/yi-462 ^ 2A ^ 2A(a2 + 862 _ 46A) '^ > 



A(A + 26-A)-2a2 _ x- + b _ 6[26(A + 26 + 3x-) - a^] , . 
2Av/X--462 2A 2A(a2 + 862 - 46A) -^ ) 

In principle, we have found the expression for T„ that we were looking for, i.e. one 
that is a function of n. However, it is somewhat unwieldy, and we will see in the next 
subsection that it can indeed be simplified. 

5. Simplification procedure 

The solution fll6p can be simplified to a compact form. To show this, let us consider 
the ordinary power series generating function of T„: 

oo 

T(z)^^T„z". (20) 

ra=0 

From ( IT6l) . we clearly have the equality 

1 — 62 1 — x+z 1 — X-Z 1 — y^z 1 — y_z 
which, after the substitution of (fTSl) and (TT9l) . becomes 

Xfz) = L±^f (22) 

Here we see the remarkable result that despite the complexity of the coefficients P± 
and Q± in (IT8l) and (iT9l) . the resulting generating function is relatively simple. We may 
proceed further and find that it may be expressed as follows: 

^ l + b^J_ f 1 1 \ . . 

^^' l-bzAz\l~x+z + b^z^ 1-X-Z + b^zy ' ^ ' 

One can recognize that the two terms enclosed by parentheses are the ordinary power 
series generating functions of the Chebyshev polynomials of the second kind [26] : 

^T , ^x sin(n + 1)6' . ,. 

Un{cos9)= \ . ^ . (24) 

We can therefore write 

-w-l^ii:h(t)-''"(lF)]-"- (-) 

n=0 

Expanding 1/(1 — bz) and using the Cauchy product gives 

-w-i^i:{i:h(t)--'(t)i^^"^" (-) 

n=0 L fc=0 
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The function T„ may thus be expressed as follows: 
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(27) 



'X- 
2b J ^""^2^ 

Given the form of the Chebyshev polynomials Uk in terms of trigonometric 
functions, we can perform the summation in the above expression over k using the 
formula for the summation of a geometric series. Following this procedure, we get 



where 



and the angles are given by 

a- ^+ 

u = arccos ■ 



cos [{n + 2)9] 
1 — cos 9 



cos [{n + 2)0] 



COS( 



arccos ■ 



X- 



26 ' " 26 

Using the half-angle formula for the cosine function, the expression for F„ 
reduces to a form 

9 



(28) 
(29) 
(30) 



Un+l ( COS - 



Un+1 I COS - 



(31) 



To make the relationship between 9 and in the above expression more illuminating, 
let us define the variable a such that 



A + 26 = 2a cos 2a. 



(32) 



Since the Hamiltonian of the system is Hermitian, its eigenvalues must be real. Hence 
A is always real. This requires a to be of the form x, ix or 7i/2 + ix, where x G M. In 
terms of a, we find that 

1 , 



where 



cos -9 
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7 = 



7e 



and cos ■ 



7e 




(33) 



(34) 



We thus have arrived at an expression for r„ in a very simple form given in terms 
of Chebyshev polynomials of the second kind: 

The relation given in fl5Sl) is the principal result of the paper. Its importance for our 
work is that we have turned the determinantal form of T„ into an expression which is an 
analytic function of n if we make use of the trigonometric definition of the Chebyshev 
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Figure 3. When atoms can only interact via next-nearest-neighbours, the chain 
decouples into to subchains, as can be seen by following the photon exchange loops 
above and below the chain. 



polynomials. If the exact roots of this equation could be found, we would have the 
explicit energy eigenvalues for single-photon states, valid for arbitrary chain lengths 
and coupling constants. Despite this complication, we can obtain the explicit energy 
eigenvalues and eigenvectors for some special cases. 

6. The case of no nearest-neighbour interaction 

We now proceed to illustrate the explicit analytical solutions of (135|) for some simplified 
cases of chains of coupled atoms. We first consider the case of no nearest-neighbour 
interaction, a ^ 0. As illustrated in Figure [3l turning off the nearest-neighbour 
interaction effectively reduces the chain into two distinct subchains that do not interact 
with each other. Let us consider a concrete example, a chain of five atoms only, as 
shown in Figure H] (a). Since the even- indexed atoms do not interact with the odd- 
indexed atoms, we are free to move them away from each other for clarity, as shown 
in Figure m (b). Notice that a photon that originated from either subchain has no way 
of being transferred to the other subchain. Thus, a five atom chain with the nearest- 
neighbour interaction ignored is equivalent to two sub-chains, one with three and the 
other with two atoms. 

This effect can also be seen in the Hamiltonian ([T]). When we put a = in the 
Hamiltonian, we obtain 

n 

H = J2^oSl + b J2 iSrSr-,2 + SrSt^,) + b J^iStSr^, + SrSt^,). (36) 

i=l jeven iodd 

It is evident that a photon in the "even" chain will never migrate to the "odd" chain, 
and vice versa. 

In practice, the procedure of removing the nearest-neighbour interaction can arise 
in two different ways: 

• For critical values of kr such that the nearest-neighbour coupling vanishes while 
the next-nearest-neighbour coupling does not. This happens for example when kr 
is equal to 4.48, 7.72, ... for the case when fi-Lvij. 
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Figure 4. (a) Next-nearest-neighbour interactions in a five atom chain, (b) The chain 
is visually separated into its two subchains; the fact that no photon lines connect them 
means that they do not interact with one another. 



• Separating all the even-indexed atoms from the odd-indexed atoms by an infinite 
distance - effectively forming two chains that are infinitely separated from each 
other. 

The other case to consider is in the limit 6-^0, which switches off the next-nearest- 
neighbour interaction. This is common practice when using the nearest-neighbour 
approximation. 



6.1. Rigorous solution in the limit of a —* 

We now illustrate the exact solution of the recurrence relation (!35|) for an arbitrary 
number of atoms with no nearest-neighbour interaction, a —>■ 0. In this limit, the 
factors 76^*" that appear in the characteristic polynomial fl35l) reduce to 

1/2 



and 



lim 7e" 



lim7e 



lim 



A + 26 + ^{X + 26)2 _ 4^2 
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v/(A + 26) /46, (37) 
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'a -f 26 - v/(A + 26)2-4a2^ 
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0. 
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We also have 

lim 2a sin 2a 
Therefore in the hmit of a - 

T 



hm v/4a2 - (A + 26)2 = i(A + 26). 







f/n+i(V(A + 26)/46)l -[[/„+i(0)] 



(39) 



(40) 



A + 26 

Equation (HOi) takes two different forms depending on whether ra + 1 is even or odd. 
If n + 1 is an odd integer, then Un+i{0) = 0, and we are left with the equation 
Un+i{\/{X + 26) /46) = 0. This equation has n + 1 roots given by a/(A + 26) /46 = 
cosk7c/(n + 2), 1 < A; < n + 1. We disregard the trivial root a/(A + 26) /46 = when 
A; = (ri + 2)/2, since the trivial root is cancelled out by the factor A + 26 present in the 
denominator of T„. The remaining n eigenvalues Ek are given by 

2k7r 



Ek = ujq — 26 cos ■ 



(41) 



n + 2' 

where 1 < A; < n + 1, excluding k = {n -\- 2)/2. Since the eigenvalues associated with 
k < {n -\- 2)/2 equal those with k > {n -\- 2)/2, we can simplify the range of k to 
1 < A; < n/2, and state that the eigenvalues are doubly degenerate. This is consistent 
with the idea that the chain should behave like two chains of length n/2 in this case. 

For n + 1 even, [f/„-|_i(0)]2 = 1, and therefore we look for the solutions of the 
equation Un+i{\/{X + 26) /46) = ±1, which are of the form 



v/(A + 26) /46 G I cos 



kn 
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1 



keIi}VJ 



kn 



cos 



n 



keh 



This implies that the eigenvalues of the one-photon states are of the form 



A G < 26 cos 



2A;7r 



n + l 



k eh\vj {2hcos 



2k7f 



n + 3 



keh 



(42) 



(43) 



where /i and I2 are sets of integers that we will determine below. 

The result (l43l) is the analytical solution for the energies of one-photon states of 
a chain of n atoms with the next-nearest-neighbour interaction only. We see that the 
energies, and thus the energy states, group into two separate manifolds. This is what 
we would expect from the above analysis of a simplified problem of five atom chain 
that in the limit of a = the atom chain effectively decouples into two separate chains. 
When the number of atoms is even, we have two sub-chains of equal length {n/2), and 
therefore doubly-degenerate energy levels. When the chain has an odd number of atoms, 
it decouples into two chains of unequal length: {n — l)/2 and (n + 1)/2. This observation 
leads us to the following restrictions on the values of k: 



A G < 26 cos ■ 



2A;7r 



n+1 



k = l,...(n-l)/2 ^U <^ 26cos 



2A; 



n 



n + 3 



k= l,...(n + l)/2 



(44) 



which gives us a total of n eigenvalues for the energies of the one-photon states. 
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6.2. Rigorous solution in the limit of b ^ 

We now illustrate the rigorous solution for the eigenvalues and eigenvectors in the other 
extreme case of no next-nearest-neighbour interaction, b ^ 0. According to (IMj) . taking 
the limit b ^ is equivalent to 7 — *> oo. 

Expanding the Chebyshev polynomials in (l35ll . we have 

T„ = ^!^ / r2n+l^n+lgi(n+l)a + . . 1 ^ _ r^n+l^n+l^-lin+l)a + . . iH _ ^5) 

2a sin 2a I '- -' '- -'J 

Clearly, the terms with the highest power of 7 will dominate the expression. We therefore 
have in the limit of large 7 

„ —i 



r22(n+l)^n+1^2(n+l) r^i2{n+l)a _ g-i2(n+l)a\l ^ j^^g-j 



2a sin 2a; 
Since 6"+^72("+i) = (a/4)"+\ we finally arrive at 

r T- ^ sin(n + l)2a 

hm Tn = a ^— . (47) 

7^00 sm 2a 

The above function clearly has n roots (since it has a Chebyshev polynomial of the 

second kind). These solutions are given by 2afc = /c7r/(n + 1) for 1 < A; < n. We thus 

obtain the energy eigenvalues 

Eu = lim (ujr, + 26 — 2a cos 2aii,) = ujn — 2a cos , (48) 

which agrees with the result of Lieb, Schultz and Mattis |lj, who made use of the 
Jordan- Wigner transform to diagonalize the Hamiltonian of the nearest-neighbour XY 
chain. An alternative derivation that does not use the Jordan- Wigner transform, but is 
restricted to the XX chain, is given in [27] . 



7. Power series solution of the general characteristic polynomial 

In this section, we will outline a power series approximation of finding the roots of the 
general characteristic polynomial (|35ll . In order to find the roots of T„, let us look 
closely into the procedure of solving the following equation 

[f/n+i (7e'")]' - [f/n+i (76-^")]' = 0. (49) 

Equivalently, this equation may be written in terms of two separate equations 

Un+l (76^") - Un+l (76-'") = 0, (50) 

and 

Un+i (7e'") + Un+i (7e-'°) = 0, (51) 

In general, the solutions a will either be real or imaginary, depending on the value 
of 7. The aim of this section is to develop a power series solution, so it would be helpful 
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to work with some property of a which does not rotate in the complex plane as 7 varies. 
Let us define 

X = ^ (arccos7e'" + arccos7e~'") , (52) 

y = ^ (arccos7e'" — arccos7e~'") . (53) 

Since we know that a must be purely real or imaginary, it is easy to see that x as 
defined above will always be real, just the property we were looking for. It is the goal 
of this section to, therefore, develop power series solutions of x for fISU]) and fl5T]) . To 
do this, it is helpful to express (150|) and (15T!) in a different form. Let us consider the 
following equation 

sm{n + 2){x + y) ^ sinjn + 2){x - y) ^ ^ 
sin(x + y) sin(x — y) 

where we obtain (ISUjl when the ± operation is replaced by a minus sign, and the converse 
for (15T|). 

Expanding out the trigonometric functions to separate x and y, we obtain 

[sin X cos y — I cos x sin y] 

X [sin(n + 2)x cos{n + 2)y + / cos(n + 2)x sin(n + 2)y] ± conj. = 0. (55) 

We have introduced the parameter / into the above expression, which is set equal to 
unity. The "conj." term is merely the first term in the expression but with the /'s set 
equal to —1. Expanding away the square brackets gives 

sin X cos y sin{n + 2)x cos(n + 2)y 

— l^ cos X sin y cos{n + 2)x sin(n + 2)y 

— / cos X sin y sin{n + 2)x cos{n + 2)y 

+ /sinxcos?/cos(r2 + 2)a;sin(ra + 2)y ± conj. = 0. (56) 

We are now presented with two expressions, corresponding to the choice made with the 
± operator. Taking the negative case causes all terms with an even power of / (or zero) 
to cancel each other out, leaving 

cos X sin y sin(n -\- 2)x cos{n + 2)y = sin x cos y cos{n + 2) sin(n + 2)y 
<^==^ tan y tan(n + 2)x = tan x ta.n{n + 2)y. (57) 

This expression is equivalent to fl50|) . We now consider the positive branch of the ± 
operator in fISB]) . In this case, the terms proportional to an odd power of I cancel each 
other out, leaving 

sin X cos y sin{n + 2)x cos{n + 2)y = cos x sin y cos{n + 2)x sin(n + 2)y 

<^=^ tan X tan(n + 2)x = tan y ta.n{n + 2)y. (58) 

We now look for power series solutions. For simplicity, we will assume that the 
chain is of even length. However, the following method will in principle work for chains 
of odd length, although the resulting power series will be more complicated. Consider 
first (jSTjl . Using the relation 

y = ^ arccos(27^ — cos 2x), (59) 
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we can express y as a function of x, and we let x to be represented by the power series 

oo 

x = Y,^kl''. (60) 

fc=0 

We next expand both sides of fl571) in powers of 7. Since we want the power series 
solution ( l60l) to be valid for all 7, we require that all the coefficients of all the powers 
of 7 vanish. Doing this gives us the following power series: 

sin Xo 2 2 cos^ a;o — 1 



x = Xq- 77-777—7 - o^: ,3,., 7 



2 cos xq 8 sin xq cos^ xq 

+ [A{n + 2)^ cos^ xo + 6 cos^ Xq - 6(n + 2)^ cos^ Xq 

—8 cos"^ xo + 2{n + 2)^ cos"^ Xq + 8 cos^ Xq — 3] 

X (48 sin^ xo cos^ xo)"S'' H , (61) 

where Xq is one of the solutions for the case 7 = 0. For n even, we have 

xo = 7T/2 + 2k7T/{n + 2), 1 < k < n/2. (62) 

A power series can be similarly obtained for (!58|) . In this case, we find that 



cos Xo 2 2 cos Xq 4 

X = Xo 7 o — 7 

2 sm Xo 8 cos xo sin xo 

+ [4{n + 2)^ cos^ Xo - 6 cos^ Xo - 10(n + 2)^ cos^ Xo 

+10 cos"^ Xo + 8{n + 2f cos^ xo - 10 cos^ xo - 2{n + 2f cos^ xo + 3] 

X (48 sin^ Xo cos^ xo)"S'' + ■ ■ ■ • (63) 

For this case, xo is simply one of the solutions to given in (!62|) . 

Given any power series, one should consider the radius of convergence. Since we 
do not have a tractable formula for [7'^]x(7) for arbitrary fc, we cannot investigate 
convergence properties. The power series presented here, however, have another 
important purpose. They provide a way to work out all the derivatives of x(7). This 
would be useful if a way could be found to formulate a differential equation to solve 
to give a closed form for x(7). The terms of the power series would then provide as 
many boundary conditions as needed at 7 = 0. It may also be possible to work out the 
generating function for the series. Some of the terms can be summed over all powers of 

7, but there are many other terms which are too intractable. Again, it is important to 
emphasize that power series terms given in this section are for chains of even length. 

8. General properties 

In this section, we shall discuss various properties of the function x(7) given by fl52|) . 
We shall not limit ourselves in this section to a power series solution. All results were 
obtained by solving ( l35l) numerically for fixed values of n. 

Working with the function x offers the convenience of investigating the energy levels 
as 7 varies from zero to infinity. The function x(7) always remains finite, as can be seen 
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in Figure O Each curve in the graph is connected with a particular energy eigenvalue 
of the system. This graph typifies the behaviour of chains with an even number of 
atoms. The graph shows which states corresponding to 7 = (a = 0, next-nearest- 
neighbour-only interaction) transition into the states associated with 7 ^ 00 (6 = 0, 
nearest-neighbour-only interaction). To show the relationship between x{'-f) and the 
energy eigenvalues, we can briefiy calculate a; (7) analytically for the cases 7 = and 
7 ^ 00. In the former case, we have (for chains of even length) that 



x(0) = h I arccos lim 7e'" + arccos lim 7e '" | 



= i arccos J(X + 2h)/Ah + \ arccos = — + - , (64) 

^ ^ 2{n + 2) 4 

where we have used the results presented in Sec. 16.11 to obtain the above result, and 
1 < A; < n/2. Therefore, we expect a total of n solution curves of x{'j) to intersect the 
7 = axis at points kir /2{n + 2) + 7r/4, 1 < A; < n/2, with two curves intersecting the 
axis at these specific points. A cursory glance at Figure [5] shows this to be the case. 

At the other extreme, we can consider the values of ^(7) in the limit 7 -^ cx). In 
this limit, we know that a ^ a^ = kTT/2{n + 1), for 1 < A; < n. We therefore have: 

lim x^j) = \ I arccos lim 7e^" + arccos lim 7e~'° ) 

y— >oo Y 7— >oo 7^00 J 



7 



1 lim (arccos 7e^"'' + arccos 7e '"'-■" 



7—* 00 



= I lim arccos (7^ — y 1 — 72e'2"'= y 1 — 'j'^e ^^"^ 

7— »oo \ 

= I arccos COS 2a;fc = a^. (65) 

Indeed, the solution curves in Figure [5] for ^(7) approach the ak in the limit of large 7; 
for convenience, the right-hand vertical axis is marked off in terms a^ = k'K/2(n + 1) to 
clearly show this. In the case of our specific graph, 7 = 2 is large enough to show this 
limiting behaviour. 

Where the lines cross, the eigenstates associated with the two levels are clearly 
degenerate. We therefore have the interesting effect that including the next-nearest 
neighbour effect can induce degeneracy, which is not present in one-photon states when 
7 — *• 00 (i.e. 6 = 0). These states are degenerate, however, when 7 = 0. When 
considering these pairs of degenerate states, we see that the highest pair has only one 
degenerate point (at 7 = 0). The second pair underneath has degeneracies for two 
values of 7, while the third pair has degeneracies occurring at three values of 7. This is 
a general trend for chains of even length. 

The solutions in terms of a also have interesting properties. These are obtained by 
simply solving fl35l) for a, which has to be done numerically for the time being, since an 
analytic solution has not been found yet for arbitrary 7. Consider Figure [6], which shows 
the absolute value of the solutions a for a 4-atom chain. We see that each solution curve 
C touches the 7 axis at only one location, 70(C). Furthermore, the curve to the left of 
7o(C) is purely imaginary; to the right, it is real. Figure [7] shows only the imaginary 
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Figure 5. The function a;(7) plotted for a 6-atom chain. 

component of a; the plot clearly supports this conclusion. This is a general trend for 
positive 7. Based on similar graphs for longer chains, these real-imaginary crossover 
points always occur for 7 < 1; we can therefore conclude that for 7 > 1, all solutions 
are real for positive 7. It is our hope that by pointing out this general feature that the 
solution curves turn through a right angle in the complex plane for some critical value 
of 7 for 7 < 1, light is shed on the possible forms an analytic solution could have. 

9. Energy eigenstates 

Given the eigenvalues Ei associated with the Hamiltonian H given by ([1]), the 
corresponding eigenvectors \tpi) for single-photon excitations can be derived as follows. 
Let the eigenvector \ipi) be given by a linear superposition 

n 

l^.) = E^^I^')' (66) 

where \j) is an eigenvector of S^ for all 1 < A; < n. The eigenvalue equation 

H\ipi) = Ei\ipi) implies that 



6(cj_2 + Cj+2) + a(cj_i + Cj+i) = E, 



j>-j, 



(67) 



for 2 < j < {n — 1) (and we have set tuo = without loss of generality for simplicity), 
and 



ac2 + fees = EiCi, 



(68) 
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Figure 6. The multi- valued function |a(7)| plotted for a 4-atom chain. There are 
four solution curves, CI, C2, C3 and C4. 



a(ci + C3) + bci = EiC2, (69) 

a(c„_2 + C„) + bCn-3 = EiCn-l, (70) 

aCn-l + bCn-2 = EiCn- (71) 

Equation fl67|) is valid for all 1 < j < n provided we impose the boundary conditions 

C_i = Co = Cn+l = Cn+2 = 0. (72) 

These boundary conditions merely force ( 167|) to reproduce ( l68l) . ( l69l) . and (ITOl) . We now 
solve the recurrence relation (1671) . Let Cj = x-'. Then ( 167|) becomes 



(73) 
(74) 

(75) 



6(x^-2 + x^'+2) + a(a;^-^ + x^+^) = E,x^ 

<^=^ 26(2 cosh^ In a: — 1) + 2a cosh In x = ii^j. 

Let z = 2 cosh In X. Then 

b{z^ -2) + az = Ei. 

This quadratic expression has the two solutions 

-a ± Va^T46(26T^ 



z± 



26 
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Figure 7. The imaginary part of a{'-f) plotted for a 4-atom chain for the four solution 
curves. 



We then have two solutions for x: 



X2 



^itarccosh (z_/2) 



We can therefore write the coefficients Cj as 



^giarccosh(2+/2) _|_ ^^-jarccosh {2+/2) 
_j_ (^giarccosh(2_/2) , j~^^-ja.iccosh{z^ /2) 



(76) 



(77) 



where A, B, C, and D are arbitrary constants whose values are derived by imposing 
the boundary conditions (!72|) . We notice, however, that we have four boundary 
conditions, thus the naive solution to this system of equations is that A, B, C, and D 
vanish. However, these four boundary conditions cannot be linearly independent; one 
of the arbitrary constants should be undetermined so that we can impose an arbitrary 
normalisation on the eigenstate. This means that we really have three independent 
equations. Without knowing the specific form of the energy eigenvalue as a function of 
a and b, however, we cannot tell which boundary condition is a trivial consequence of 
the other three. Therefore, we stop the derivation at this point. We have laid out most 
of the calculation so that once the energy eigenvalues can be completely determined, 
the remaining steps to determine the eigenstates are simple. 
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10. Conclusion 

We have derived tlie cliaracteristic polynomial of the eigenvalue equation for a chain 
of atoms with up to the next nearest-neighbour coupling. The exact form of the 
polynomial is given in terms of the Chebyshev polynomials of the second kind and 
is valid for an arbitrary number of atoms and coupling strengths. We have illustrated 
the exact analytical solutions of the characteristic polynomial for two simplified cases 
of no nearest-neighbour coupling and no next-nearest-neighbour coupling. 

We have also derived the power series solution of the general characteristic 
polynomial. We have shown that there are patterns as to how many times a pair 
of eigenstates become degenerate. We have also outlined a simple approach by which 
the energy eigenstates for single excitations may be found. 

We hope that our presentation provides tools for further advancement towards an 
analytic solution for the energy eigenvalues of arbitrary excited linear chains and in the 
analytic studies and understanding of the energy structure of linear chains of coupled 
atoms and also to simplification of the mathematics involved. This alone hopefully is 
advantageous since it helps reduce the number of operations a computer has to perform 
to find the energy eigenvalues. 
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